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ABSTRACT 


Pressure and velocity data were collected in a full scale jet 
engine test cell in order to validate the predictive accuracy of a 
two-dimensional and axisymmetric primitive variable computer code. It 
was found that the model reasonably predicted the velocity profiles in 
the augmentor tube. Inaccuracy increased at higher engine thrust settings 
at positions far downstream in the augmentor tube. Predicted pressure 
profiles were reasonable but the magnitudes were in considerable error 


at high flow rates. 
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I. INTRODUCTION 


Jet engine test cells are in wide use today in military and com- 
mercial aviation maintenance programs. Their purpose is to provide 
for testing of jet engines under conditions that approximate the in- 
stalled environment, thus enabling the responsible authority to veri- 
fy the engine's capabilities throughout its operational envelope 
prior to returning it to service. However, due to exhaust pollution 
problems inherent to their operation, the use of these test facilities 
has recently encountered a considerable amount of criticism and even 
judicial action in the State of California. 

The necessity for testing an engine after overhaul is apparent. 
However, the pollution problems that result when testing today's high 
power, high mass flow engines in present test facilities must also be 
considered. 

"The typical test cell incorporates an inlet, a horizontal test 
section and vertical exhaust stack. The engine to be tested is 
normally mounted near the center of the cell to allow the development 
of a nearly uniform engine inlet velocity profile. The engine ex- 
nausts into an augmentor tube which entrains additional air for ex- 
naust gas cooling and dilution. The quantity of this secondary air 
is crucial to proper engine testing and test cell performance" 


[Ref. 1]. The ratio of entrained air mass flow rate to engine mass 
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flow rate is known as “augmentation ratio’. An augmentation ratio that 
is too high could result in inaccurate engine performance measurements 
and in exceeding the test cell structural limits due to pressure 
gradients and excessive cell pressure reduction. Insufficient secondary 
air could allow engine exhaust gas recirculation to the engine inlet and 
cause excessive opacity of visible emittants. It is therefore necessary 
to develop mathematical models which can be used to predict the flow 
fields within turbojet test cells. These models could then be used as 

a cost effective means for determination of optimum test cell designs 
for meeting the military engine testing needs while minimizing noise and 
chemical pollution effects on the local environment. 

With this in mind, considerable experimental and computer modelling 
work has been done at the Naval Postgraduate School to further the 
understanding of the flow field and detailed operating characteristics 
encountered in typical test cells. This knowledge is currently needed 
to facilitate the necessary modifications to present test cells and for 
the design of new cells to accommodate the high technology engines of 
the future. 

The two-dimensional modelling of the flow within turbojet test 
cells was initiated at NPS by Hayes and Netzer [Ref. 2]. They adapted 
the two-dimensional, w-w model of Spalding, et al [Refs. 3, 4] to the 
test cell geometry and operating conditions. This work was continued 
by Speakman, Hayes and Netzer [Ref. 5]. The model made use of the 
fundamental equations of conservation of mass, momentum and energy. 


The dependent variables were stream functions (wv), vorticity (w), 





turbulence kinetic energy (K), the turbulence dissipation rate (ce), and 
the temperature (T). As a consequence of the w-w formation, pressure 
(and to a large extent velocity) is removed from the governing equa- 
tions. These second order elliptic partial differential equations 

were reduced to "finite-difference," non-linear algebraic equations and 
solved simultaneously with an iterative procedure. Results of this 
computer simulation appeared to be quite accurate at low subsonic en- 
gine exhaust Mach numbers (M < 0.6), but pressure predictions were un- 
reliable at high subsonic or sonic exhaust velocities. These compari- 
sons between predictions and experiment were made by Walters and Netzer 
[Ref. 6 | using a one eighth scale model of an NAS Alameda test cell. 
Other problems related to the ~-w formulation were: 1) restriction to 
constant density flows, or to flows in which density varied only with 
temperature; 2) boundary conditions were difficult to specify; 3) con- 
Siderable difficulty was experienced in obtaining converged solutions, 
especially for non-uniformly spaced grids and high flow rates, and 4) 
the ~-w model is not easily extended to three dimensional flows. 

To help alleviate the above modelling difficulties Stevenson and 
Netzer [Ref. ip utilized a primitive variable (pressure-velocity) model. 
This computer model (adapted from the CHAMPION 2/E/FIX computer program 
developed by Pun and Spalding [ Ref. 7| extended the modelling capabili- 
ties to more complex geometries and to high subsonic and sonic engine 
exhaust velocities. It also is readily applicable to variable density 
flows. The earlier results of the w-w model and those of the u-v-p 
model were compared to the empirical data obtained with the subscale 


test cell [ Ref. 6 |. Good agreement between the two methods was found 
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for low subsonic engine exhaust conditions. The primitive variable 
computer model produced reasonable results for the flow field downstream 
of the engine exhaust up to the sonic exhaust condition. 

The primary objective of this investigation was to adapt Spalding's 
CHAMPION 2/E/FIX primitive variable computer model (as modified by 
Stevenson) to the geometry and flow conditions of an actual operational 
full scale turbojet test cell and to test that model with empirical data 
from the same cell. The turbojet test cell chosen was one located at 
NAS Alameda, operated by NARF Alameda. The engine being operated at 
the time of data acquisition was the General Electric TF34-GE-100, cur- 
rently used on the USAF A-10 attack aircraft. The major computer model 
modifications required were: 1) to increase the geometry to the full 
scale, and 2) to change the uniform engine exhaust velocity profile to 
a non-uniform profile which adequately represented the primary and 


secondary flows from the TF-34 turbofan engine. 





IIT. EXPERIMENTAL INVESTIGATION 


A. INTRODUCTION 

The turbojet test cell chosen for this investigation was test cell 
No. 13 located in Building 372 at NAS Alameda. It is operated by NARF 
Alameda as the prime means of evaluating the USAF TF34-GE-100 turbofan 
engine following required normal periodic maintenance. 

The TF34-GE-100 is a turbofan engine in current use on the USAF 
twin engine A-10 attack airplane. The engine is 2.54m long and produces 
a maximum rated thrust of 9065 Ibf [Ref. 8]. The bypass ratio is 6.22:1 
with maximum rated airflow at 151.0 kg/sec. The exhaust nozzle is a 
fixed area, convergent type with the cold stream being 0.474 sq. m. in 
area and the hot stream 0.148 sq. m. in area. While operating at maxi- 
mum power the core exhaust gas is estimated to have a velocity of 305 m/ 
sec and a temperature of 393°C [ Ref. 8]. 

To provide the data for model validation it was decided to measure 
velocity profiles within the test cell and at two locations within the 
augmentor tube as well as the axial pressure distribution throughout 
the test cell. To most accurately determine velocity profiles, mul- 
tiple pitot-static tubes with associated stagnation temperature measure- 
ments would have been required. However, due to cost and time restraints 
involved with the manufacturing and installation of such equipment, a 
compromise approach to this problem was taken. Two stagnation tube 
rakes, each fitted with two thermocouples, were installed in the aug- 
mentor tube. Also, thirty two static pressure ports were installed in 


the augmentor tube wall. Velocity measurements were also made within 
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the test cell using hot wire velocimeters. Details of the data acquisi- 
tion as well as the method of data reduction are presented below. 
peeeees! CELL 

The turbojet test cell utilized in this investigation is typical of 
test cells currently used in military and commercial aviation maintenance 
programs. It was a concrete blockhouse type structure 43m long, 9m high, 
and 9m wide (Fig. 1). It had a horizontal inlet, which was retangular in 
cross section, and exhausted into a vertical stack. Not shown in Figure 
] are the noise suppression devices which were located over the cross 
Sections of both the test cell inlet and back wall. Located 15m from 
the cell inlet was the test platform which stood 3m above the floor of 
the cell. The jet engine test stand and the first section of the aug- 
mentor tube were located on this platform. The TF-34 was mounted on the 
test stand with an attached, fixed geometry bellmouth type inlet, along 
with necessary start-up and test equipment to monitor engine performance 
in the cell control room. This test stand was movable, so that the 
distance between the exit plane of the engine and entrance plane of the 
augmentor tube could be varied. During this investigation this distance 
was fixed at 1./m. 

The first section of the augmentor tube was located on the test 
platform. As shown in Figure 1 it too was movable. It was a steel 
pipe with 9.5mm thick walls and was approximately 1.8m in diameter and 
1.6m long. The next section was also a constant diameter pipe, 5.2m 
long. During this investigation, the pipe protruded only 5m beyond 
the first section. The radial gap between these two sections was, on 


the average, l6mm. The next section was 3.5m long and diverged in 





diameter from approximately 1.8m to 2./m. At the downstream end of the 
latter section was a water-injection ring. Following a 0.3lm gap was 
the next section: a constant 3.05m diameter pipe, 13.1m long. In 
this section were three staggered rows of four noise suppression de- 
vices. Each unit was a 3.2m long cylinder, 0.6m in diameter. 
C. INSTRUMENTATION AND DATA ACQUISITION 

Two pitot tube rakes were utilized in this investigation. Each 
rake was constructed from 102mm wide, 9.5mm thick, 302 stainless steel. 
The cross-like configuration used is depicted in Figs. 2 and 3. The 
first rake was located 127mm from the entrance plane of the augmentor 
tube. The tips of the pitot tubes being 25.4mm from the entrance. 
The vertical member of this rake was |.67m high and had seven pitot 
tubes built into it. Each horizontal member was 0.822m long with three 
pitot tubes inserted. Horizontal and vertical members were attached as 
Shown in Fig. 4. To hold the rakes in position, 50.8mm high tabs were 
welded into the augmentor tube at the 12, 3, 6, and 9 o'clock positions. 
The rakes were then attached to the tabs as shown in Fig. 5. The down- 
Stream rake was constructed exactly the same as the first except that 
the vertical length was 1.70m and each horizontal member was 0.833m 
long. This second rake was installed 5.05m downstream from the aug- 
mentor entrance plane (127mm upstream from the start of the diverging 
section of the tube). The 6.35mm diameter stainless steel tubing used 
for the pitot tubes had a wall thickness of 1.57mm. Each tube pro- 
truded 101.6mm forward from the rake. The tubes were silver-soldered 
on the back side of the rake and routed from the augmentor tube at the 


aft end of each attachment tab through holes drilled in the augmentor. 
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Tygon tubing (3/16 in. ID x 1/16 in. wall) was then routed from pitot 
tube to a Scanivalve located in a standing cabinet next to the jet- 
engine test stand. The output of the Scanivalve was routed through the 
floor of the cabinet and underneath the test platform to the recording 
apparatus in the control room (Fig. 6). In addition to the pitot tubes, 
thermocouples were mounted on each rake as shown in Fig. 2. The thermo- 
couple wire was routed from the augmentor with the pitot tubing, then to 
ice bath references. The thermocouple data were obtained from a digital 
voltmeter. 

Static pressure ports were located at the 2, 5, 7, and 10 o'clock 
Meese vons at axial locations 0.152m, 0.305m, 0.610m, 1.212m, 3.96m, 
5.18m, 6.10m, and 7.32m downstream of the augmentor inlet. Circumferen- 
tial ports were used to determine whether or not the flow was axisym- 
metric. Tygon tubing was attached to Swagelock fittings at each port 
and routed to a second Scanivalve in the standing cabinet. 

In order to determine the velocity of the secondary (test cell) air, 
velocimeters were located 0.6m and 1.8m above the TF34 engine. The 
velocimeters used were Datametrics Air Flow Meters, commonly used to 
measure the velocity in air conditioning ducts. Prior to their use, 
calibration was performed at the NPS using a small wind tunnel. The 
resulting calibration curve is shown in Fig. 7. The output of these 
velocimeters were routed underneath the jet engine test stand to the 
Standing cabinet and then to the control room as described above. In 
the control room the output signals were connected to a digital volt- 
meter. 

The instrumentation apparatus set up in the control room (Fig. 6) 


consisted of a locally manufactured Scanivalve controller, two Doric 





digial voltmeters and two Honeywell “Electronic 196", 2 pen strip chart 
recorders. 
D. EXPERIMENTAL PROCEDURE 
Prior to actual data acquisition the following steps were performed: 
1) The 30.5 cm radial gap between the second and third augmentor 
Sections was blocked off with 1.9 cm thick plywood. This was 
done to prevent any secondary cell air from entering the aug- 
mentor downstream of the data acquisition equipment. It should 
be noted that an airtight seal was not attempted. 

2) All static and dynamic pressure ports were pressure checked to 

insure proper connections and identification. 

3) All thermocouples were checked for proper rake attachment, 

wiring and identification. 

The actual data acquisition process began before anticipated due to 
inadequate thermocouple installation. After initial TF34 start up and 
a brief run time (during which the recording instruments were calibrated) 
the test cell operators discovered an oi] leak in the TF-34 and conse- 
quently shut it down for repair. This provided an opportunity to 
inspect the test equipment. At this time it was discovered that the 
thermocouples were not suitably constructed to withstand a complete 
series of tests. Therefore, it was decided to record termperature data 
during the initial TF-34 bearing breakin tests at various power settings 
for as long as the thermocouples remained in place. Reliable tempera- 
ture data were only obtained for maximum, 90%, and idle power settings. 


The data appears in Table I. 





POWER 


CORE FAN 


SETTING RPM RPM 


100% 


90% 


IDLE 


TABLE 


LUeZe leletsy// 


16120 5803 


Prs50 0 685 


THRUST 
LBS 


8611] 


Soa7 


382 


RAKE 


FORWARD 


AFT 


FORWARD 


AFT 


FORWARD 


AFT 


THERMOCOUPLE 


CENTER 
LOWER 
CENTER 
LOWER 


CENTER 
LOWER 
CENTER 
LOWER 


CENTER 
LOWER 
CENTER 
LOWER 


TEMPERATURE 


46] 


156 


498 


105 


I: Thermocouple Data Obtained During Initial Engine Testing 
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After the brief bearing breakin tests the TF-34 was put through a 


series of five minute duration tests at a variety of power settings. 


When a particular power setting was reached two minutes were allowed 


to insure steady-state operation and then the following data were re- 


corded: 


1) From the test cell instrumentation: 


a 


at ‘-) Oo 


( 


) 
) 


i 


ambient temperature 
ambient pressure 
cell pressure 

TF-34 core RPM 
TF-34 fan RPM 

TF-34 thrust 


2) From NPS installed equipment and instrumentation: 


a) 
b) 
e} 


Static pressures along the augmentor wall 
pitot tube pressures at each rake 


temperatures at each rake 


Each 48 channel Scanivalve was set up so that port number one was 


Open to the cell pressure. A fixed pressure (set using a mercury 


manometer) was input to port number two in order to provide a refer- 


ence calibration each time a cycle was begun. During each run one 


Scanivalve was cycled using the Scanivalve controller, data recorded 


On the strip chart, connections switched, and the same procedure 


followed for the remaining Scanivalve. The pitot rakes were labeled 


as shown in Fig. 8, and static ports as in Fig. 9. The recorded data 


are presented in Tables II through IV, the letter "A" referring to the 


forward rake, "B" to the aft rake, and letters "C through J" to the 
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static port stations from front to rear. The data obtained from the 
intermediate power setting run was not utilized due to lack of reliable 
temperature information at that thrust level. 

There were three specific problems encountered during the data ac- 

quisition process: 

1) Failure of thermocouple attachments. This resulted in lack of 
temperature data during each run that pressure data were ob- 
tained. However, this was overcome by correlating (by thrust 
setting and fan and core speeds) the temperature data recorded 
at the start of the TF-34 test series with the pressure data 
obtained later. 

2) During the entire test series a loud resonant frequency was 
heard. This initially forced engine shut down to allow for 
cell inspection. No damage was discovered and all instrumen- 
tation appeared normal so testing was continued as planned. 
However, upon removal of test equipment it was discovered that 
two of the welds holding the aft rake in place had broken. It 
should also be noted that at the same time it was discovered 
that two pitot tubes had detached from the same rake. The 
latter apparently occurred during the last test since no 
unusual change occurred in the recorded data. 

3) Due to electrical connection difficulties the velocimeters were 
not installed during the engine testing periods when other data 


were obtained. When this problem was corrected (after the pres- 


sure data were obtained) it was discovered that the wooden blocks 


2] 





RUN NO.: 1 POWER SETTING: 100% TIME: 1305 
PITOT RAKE (p~Poey, ) STATIC PORT (p-Poc,, ) 
SCANIVALVE ie) SCANIVALVE RES 
cuanne, PORT I.D. (N/m“x10°~) CHANNEL PORT 1.D. — (N/m°x10 ) 
] Zero ] Zero 
2 REF PRES 2 REF PRES 
3 Al 34. 30 3 CI -~14.196 
4 A2 39 .93 4 GZ -13.776 
5 A3 3.56 5 C3 -13.803 
6 A4 -0.63 6 C4 -14.266 
7 A5 41.11 7 DI = ena 
8 AG 5.97 8 D2 = 2a 
9 AT -0.48 9 D3 -12.956 
10 A8 39.58 10 D4 -13.190 
1] AQ 3.70 11 Ey -9.694 
12 A10 =084 12 E2 = 9370 
13 All 40.19 13 F3 ~9 432 
14 Al2 2.68 14 E4 -9.860 
15 A13 -0.48 15 Fy -8.226 
16 Bl 26.84 16 F2 -8.226 
7 B2 19.81 17 F3 MY YAS 
18 B3 8.39 18 F4 -8226 
19 B4 3.30 19 Gl -6 433 
20 B5 23.24 20 G2 -6.661 
21 B6 10.36 2] G3 -6 .530 
22 B7 3.12 22 G4 -6 433 
23 B8 18.71 23 H1 -6.268 
24 B9 7.51 24 H2 -6.171 
25 B10 1.28 25 H3 -5.875 
26 Bl] 15.59 26 H4 -6.33 
27 B12 5.76 27 1] ~4.144 
28 B13 0.04 28 12 -4.013 
29 13 ~3.820 
30 14 rai (Ory) 
31 om le 4] 
32 JZ -1.958 
33 Je “2527 
34 J4 my ie7, 
FAN CORE 
RPM: 6701 17360 
THRUST: 8685 
(LBS) 
CELL RAKE A RAKE B CELL AMBIENT 
TEMP: 17.8 PRESSURE: 1.0129x10° 1.0141x10 
ec) (N/m) 
TABLE II: EXPERIMENTAL DATA, 100% POWER 
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RUN NO: 3 


PITOT RAKE ( 


SCANIVALVE 
CHANNEL 


FAN 

RPM: 5933 
THRUST: 
(LBS) 


PEEL 
lhe 8 


TEMP: 
-C) 


6304 


POWER SETTING: 


P-Pcet? 

PORT 1.D.  (N/m*x107 
Zero 

REF PRES 
Al 22.64 
A2 30.86 
A3 0.49 
Ad -0.39 
AS 29.04 
AG 4.18 
AT -0.19 
Ag 30.60 
AQ 2.58 
A10 =), 2 
All Ay 7 
Al2 3.65 
A13 0). 
BI 20.64 
B2 16.86 
B3 8.18 
B4 1.60 
BS 17.30 
B6 8.4] 
B7 3.16 
Bg 13.96 
B9 5.83 
B10 1.38 
B11 [27 
B12 3.96 
B13 0.98 

CORE 

16403 

RAKE A RAKE B 
TABLE III: 


90% TIME: 1327 
STATIC PORT (p-Poc)) ) 

3) aes latte PORT I.D.  (N/m*x107°) 
1 Zero 
2 REF PRES 
3 C1 9,29 
4 C2 -10.08 
: C3 EGns6 
6 cA 9.29 
7 D1 9.23 
8 D2 -9.23 
9 D3 By 
10 D4 9.45 
1 F -746 
12 E2 -6.70 
13 £3 BF 
14 Eq Gs 
15 Fy i 
16 F2 5.62 
7 F3 Be oe 
18 FA nS Be 
19 GI Pia 
20 G2 -4.41 
21 G3 4.4) 
22 GA -4.26 
23 M1 -4_08 
24 H2 -3.87 
25 H3 Ba 5 
26 Hd 4.12 
27 1 -2.76 
28 12 2.79 
29 13 2.45 
30 14 -2.61 
31 J -1.10 
32 32 -1.40 
33 J3 -1.40 
34 34 Sed 

CELL , AMBIENT 
PRESSURE: 1.0136x10° 1.0144x10 


(N/m*) 


EXPERIMENTAL DATA, 90% POWER 
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RUN NO: 4 POWER SETTING: IDLE Woges 
PITOT RAKE (P-Poe,) SPA PORT (D-Peey | ) 
SCANIVALVE 2 -3 SCANIVALVE 2 
CHANNEL PORT I.D. (N/m™x10 ~) CHANNEL PORT I.D. (N/m 
] Zero 1 Zero 
2 REF PRES Z REF PRES 
3 Al 1.44 3 Cl -0 
4 A2 Ze56 4 EZ -0 
5 A3 0.45 5 oe: -0 
6 A4 acy 6 C4 -Q 
7 A5 252 i, DI -0 
8 A6 0.4] 8 D2 -0 
3 A7 Ot 9 D3 -0 
10 A& a N6| 10 D4 -0 
1] A9 On27 1] El -0 
12 Al0 0.14 FZ EZ -0 
18 All 2.48 Vs ES -0 
14 Al2 0.23 14 E4 ~Q 
115 Al3 0.14 15 F] -0 
16 Bl ilieos 16 F2 -0 
Ay B2 1.48 7 FS -0 
18 B3 0.54 18 F4 -0 
19 B4 0.50 19 GI -0 
20 B5 O22) 20 G2 -0 
21 B6 1.35 21 G3 -0 
22 B7 35 22 G4 -0 
23 B8 0.36 23 H1 -Q 
24 B9 Pele 24 H2 -0 
25 B10 0.54 25 H3 -0 
26 Bl] 0.32 26 H4 -0 
27 B12 ls IZ 2] te -0 
28 B13 0.41 28 {2 -0 
29 [3 0 
30 14 0 
oH J] -Q 
Sic J2 0 
os J3 0 
34 J4 0 
FAN CORE 
RPM: 1801 11404 
m7RWST: 410 
(LBS) 
CELL RAKE A RAKE B CBE 5 AMBIENT 5 
HEMP: 18.3 PRESSURE: 1.0144x10 1.0144x10 
ec) (N/m?) 
TABLE IV: EXPERIMENTAL DATA, IDLE POWER 
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| Sa 


x107°) 


.64 
di 
74 
va 
Wl 
74 
74 
wl 
nels, 
49 
ore 
5218 
40 
43 
40 
43 
as 
BA) 
s20 
SAS 
.28 
.18 
.18 
les 
.03 
503 


m0) 





at the intersection of the second and third augmentor sections 


had partially detached. This caused the velocity of the cel] 
air to probably be somewhat higher during the tests using the 
velocimeters than it was during the pressure/acquisition runs 
when the blocks were in place. The velocimeter data are shown 
in Table V. These data were used only to obtain the general 
shape of the velocity profile above the engine. 
E. DATA REDUCTION 
The a data that appear in Tables II through IV came directly 
from the strip chart recordings, or, in the case where the reading was 
too small, from the digital voltmeter. The calibration scale on the 
strip chart (psi/division) was obtained from a 1.05 in. Hg reference in- 
put to Channel No. 2 of both Scanivalves. Using this information, 
axial pressure profiles along the four circumferential positions were 
plotted and an average pressure profile was then estimated for each 
power setting (Figs. 10, 11, and 12). These profiles were then used to 
obtain the static pressure at the aft rake location (it was assumed 
that pressure was uniform across the augmentor at the aft rake). 
The measured stagnation pressure at each pitot tube together with 
the static pressure at the rake were used to obtain a Mach number pro- 


file using the expression 
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Velocity RPM Thrust 


(m/sec) (lbs) 
Upper Lower 
Velocimeter Velocimeter Fan Core 
ee 5 | oe 1789 T2265 520 
Z03 3.05 6780 17285 9072 


TABLE V: VELOCIMETER DATA 
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with y assumed equal to 1.4. The thermocouple readings were plotted and 
used to estimate an assumed axisymmetric stagnation temperature profile 
for the aft rake at each power setting (Fig. 13). The estimated stagna- 
tion temperature profile and local Mach number were then used to deter- 
mine a static temperature for each pitot tube location. Using this in- 
formation and the properties of air at standard temperature and pressure, 
accoustic speed and velocity were found for each pitot tube position. 
Static temperature for the aft rake is plotted as a function of radius 

in Fig. 14, velocity profiles are shown in Figs. 15, 16, and 17. 

The static pressure of the front rake was estimated by extrapolating 
the average pressure profile to the augmentor inlet position. This was 
required due to the fact that the first set of static pressure ports 
was installed 15.2 cm aft of the augmentor entrance plane. This static 
pressure was assumed to be uniform across the augmentor for each power 
setting. Mach number profiles were then obtained as above. At the 
forward rake y was not assumed to be uniform across the augmentor. In- 
Stead, y was taken to be 1.33 at the centerline position and at the four 
pitot tubes located circumferentially 25.4 cm from the center. At the 
remaining positions y was assumed to be 1.4. This property change was 
made since heat induced discoloration on the forward rake was observed 
to occur from the core jet to a radius of approximately 30 cm. Stagna- 
tion temperature profiles were estimated using the measured temperatures 
at the two rake locations and are shown in Fig. 18. Using this stagna- 


tion temperature distribution with the Mach number profiles, accoustic 


Ci 





Speed and velocity were determined as above. The resulting velocity 


profiles for each power setting are shown in Figs. 19 through 21. 
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eco UNER MODEL 


A. INTRODUCTION 

The computer program used to model the turbojet test cell of this 
investigation was the CHAMPION/2/E/FIX computer program developed by 
Pun and Spalding [Ref 7] was modified by Stevenson [Ref. il. This is a 
primitive variable (pressure-velocity) model for two-dimensional, turbu- 
lent recirculating flows. The equations are cast into finite difference 
form and solved using a line-by-line iteration method. A detailed dis- 
cussion of the model can be found in Refs. 1 and 7. 
pee PHYSICAL MODEL 

The major modifications made to the computer program as utilized by 
Stevenson [Ref. i| were to change the geometry of the test cell and to 
incorporate the engine exhaust velocity profile of the TF-34 turbofan 
engine. 

The turbojet test cell, TF-34 engine and augmentor dimensions used in 
the computer model are shown in Fig. 22. There were several differences 
between the configuration used in the computer model and the actual test 
cell as shown in Fig. 22: 

1. The computer model was axisymmetric and therefore required the 

engine to be located on the axis of symmetry of the test cell. 

In the actual cell the engine is approximately 0.61m closer to 
the deck than to the overhead. Thus, it can be expected that the 
predicted palo distribution of the secondary cell air would 


be somewhat different than the actual distribution. In addition, 


a 





the test stand and supporting deck probably distorted the velocity 
profile at the augmentor inlet. 

The exhaust of the TF-34 turbofan engine was modeled as two con- 
centric cylinders (gas generator core and fan exhausts). In 
reality, each exit incorporates a converging nozzle. For simpli- 
city in the two-dimensional model, the actual exit diameter of 
each duct was utilized over the entire engine length (Figs. 23 
and 24). 

In the computer model the aft cell wall was taken to coincide 
with the augmentor inlet. This wall was actually approximately 
llm farther back. The effect of this modification was to elimi- 
nate the recirculation zone above the augmentor tube. However, 
the effects of this low velocity recirculation zone above the 
augmentor flow field have been shown to be insignificant [ Ref. 5 |. 
The computer model utilized a constant diameter augmentor tube 
20.5m long instead of the varying area tube construction shown 

in Fig. 1. The second rake was not positioned farther down- 
Stream to avoid complications arising from the divergent 
character of the second augmentor section and the flow distur- 
bance expected to result from the water injection ring and 

noise suppression devices. The augmentor tube length was 
modelled as 20.5m versus the actual 5.05m so that the boundary 
conditions at the exit plane would not influence the calcula- 
tions at the position of the second rake. 

The test cell was square in cross-section. This was modelled as 


a concentric cylinder with the same cross-sectional area. 
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The grid spacing for the finite difference numerical solution of this 
problem was crucial. A 30 by 30 grid was used originally without obtain- 
ing convergence. Consequently the grid was expanded to a 40 by 30 to 
better accommodate the much greater length of the axial direction. The 
recommendation by Gosman et al | Ref. 3 | that Successive spacing should 
not increase by more than a factor of approximately 1.5 was adhered to. 
C. DATA INPUT 

Using the empirical velocity profiles from Figs. 15 through 17, an 
average velocity for each radial position at the aft rake was determined. 
These average velocities and the centerline velocity were then plotted 
(Fig. 25) and graphically integrated to obtain a one-dimensional bulk 
velocity at the aft rake for each power setting. An identical procedure 
was followed to obtain an average static temperature (Fig. 14) at the aft 
rake. These results were then used (together with the standard properties 
of air and known augmentor tube cross sectional area) to calculate the 
total mass flow rate through the augmentor tube. 

Using the fan and core RPM, thrust, and engine inlet temperatures 
that appear in Tables II through IV, fan and core mass flow rates for 
the TF-34 were obtained from the General Electric Corp., Lynn, Mass. 

It should be noted that the flow rates provided by General Electric 
were for "an ideal" engine and may not exactly apply to the specific 
engine used in this investigation. 

With the above flow rates it was then possible to calculate the 
mass flow rate through the test cell. Then, using the temperature and 
pressure in the test cell (and the known cell cross sectional area) the 


average one-dimensional velocity of the cell air was determined. 
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In this investigation no instrumentation was available for measuring 
Static pressures and temperatures at the TF-34 core and fan tailpipe ex- 
hausts. Turbine exit temperature as a function of thrust is presented in 
Ref. 8. 

This temperature was used as the tailpipe stagnation temperature in 
order to calculate nozzle exhaust temperature. For the idle power setting 
the Mach number at the core exit was assumed to be approximately 0.5 and 
y to be 1.33. Using a turbine exhaust temperature of 1560°R resulted in 
a nozzle discharge static temperature of 1498°R. For the 90% and full 
power cases the nozzie flow was choked and the turbine exhaust tempera- 
tures were 1694°R and 1922°R respectively. This resulted in exhaust 
temperatures of 1453°R and 1649°R respectively. Fan exhaust temperatures 
were available from the temperature profiles of Ref. 8. Actual cross 
Sectional areas of both exhaust nozzles were measured on the engine. 
Static pressure at the exhaust nozzles was taken to be the static pres- 
Sure Of the cell in all cases. This static pressure estimate is 
accurate for the fan exhaust in all cases and the core exhaust at the 
idle power setting. However, the static pressure of the core exhaust 
for the maximum and 90% power settings was actually the critical pres- 
Sure. 

Using the above information the one-dimensional continuity equation 
was used to calculate the bulk velocities for core and fan exhausts at 
all power settings. The computer input data as well as the CPU time 
required to obtain a converged solution for each power setting are shown 


in Table VI. 
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POWER SETTING TEMPERATURE VELOCITY CELL-PRESSURE CPU TIME 
( ) 


(°K) (m/sec) N/m™-x10- IBM 360/67) 
Pee 29 OFeZ25 TOS Tks, 
IDLE FAN 344 96.04 
CORE 832 S38) 196: 
GEE 05) Zh 101.37 79 
90% FAN 305 174.2 
GORE 807 ZO8eo 
CEE 290 Bee Gllese 150 
1007 FAN 306 202.9 
CORE 916 354.1 


TABLE VI. COMPUTER INPUT DATA 
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Ri  UPSSUSSTON OF RESULTS 


A. INTRODUCTION 

The results of this investigation are presented in the following 
manner: 

1) General characteristics and accuracy of the experimental data 
and sensitivity of the calculated experimental velocity pro- 
files to assumptions made for the unknown/unmeasured experi- 
mental parameters. 

2) A comparison of experimental data to computer model predictions 
and the sensitivity of the predictions to uncertainties in the 
data. 

B. EXPERIMENTAL DATA 

Experimental data are presented in Tables II through IV. The re- 
sulting calculated velocity profiles for the fore and aft rakes at each 
power setting are presented in Figs. 19 through 21 and Figs. 15 through 
17, respectively. Figs. 10 through 12 show the static pressure profiles 
along the augmentor wall for each power setting. 

The pressure profiles show only a small variation with circumferen- 
tial location. In all cases the minimum pressure occurred very close 
to the inlet plane of the augmentor, whereas the maximum pressure in 
the instrumented portion of the-augmentor at no time exceeded atmos- 
pheric. As expected, the minimum pressure at the augmentor inlet de- 


3 N/m° 


creased as the engine flow rate increased, going from -0.8x10> 
at the idle power setting to -14.0x107? N/m? at the maximum power condi- 
leon. 
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The velocity profiles at the aft rake indicate that the flow was not 
axisymmetric in the augmentor tube. The most asymmetric condition oc- 
curred at the idle power setting. The idle profile was also much more 
fully developed than the others, indicating that the flow at idle was 
the only one that entered the diverging section of the augmentor tube in 
a reasonably well mixed condition. There also appeared to be a general 
trend toward higher mass flow in the quadrant defined by the 12 and 3 
Seelock positions. 

In all cases the velocity profiles at the forward rake showed the 
expected high velocity near the center with the lower fan and augmentor 
air velocities toward the wall. It should also be noted that the spread 
in data at the forward rake was much less than that at the aft rake. 

In all cases the engine-exit-to-augmentor-inlet spacing was 1.7m. At 
the idle power setting the augmentation ratio was 2.39:1. This ratio 
dropped to 2.22:1 for the 90% power case and rose to 2.42:1 for the maxi- 
mum power setting. 

In summary it can be said that the static pressure remained quite 
uniform across the tube, that the flow at the aft rake was not axisymmet- 
ric, and that at the front pitot rake the flow was quite symmetrical. 

For the 90% power setting the stagnation pressure data at the aft rake 
was estimated to have an uncertainty of 0.17% which, in turn, resulted in 
a 0.25% uncertainty in the calculated velocities. Similarily, the sta- 
tic pressure readings were examined and found to have an uncertainty of 
approximately 2% near the aft rake. This uncertainty in static pressure 
resulted in approximately a 5% uncertainty in the velocities at the aft 


rake. 
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The possibility that the total temperature profile at the aft rake 
was not that estimated in Figure 13 was also examined. A somewhat dif- 
ferent, but plausible, total temperature profile was plotted and is shown 
in Fig. 26 for 90% power. This change in the Ty profile resulted in a 
4.7% increase in the bulk static temperature and a 1.3% increase in the 
bulk static temperature and a 1.3% increase in the bulk velocity. These 
changes would result in a 4.8% decrease in the test cell mass flow rate 
(assuming engine flow rates were known). 

Finally, the assumption that the specific heat ratio was not 1.4 at 
the aft rake was examined. Changing y to 1.35 resulted in insignificant 
changes in the cell flow rate. 

The experimental data for the forward rake were also examined to dis- 
cover how uncertainties in static pressure and stagnation temperature 
profiles would affect the predicted velocities. Possible errors in the 
extrapolated pressure at the forward end of the pitot tubes resulted in 
an uncertainty of only 4.5% (Fig. 27) in the velocities. The velocity 
profile is sensitive to the assumed stagnation profile in much the same 
manner as discussed above for the aft rake. 

C. COMPUTER MODEL PREDICTIONS 

The pressure and velocity profiles that resulted from computer 
modelling are presented with the actual experimental data for each 
power setting in Figs. 10 through 12, 15 through 17, and 19 through 21. 

Sensitivity of the computer model predictions to specified test 
cell mass flow rate and engine mass flow rate were examined. A TO% de- 
crease in engine mass flow rate was found to alter the computer model 


results an insignificant amount. This is attributable to the lesser 
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contribution of the engine to the total flow rate through the augmentor. 
A 10% decrease in the specified test cell mass flow rate decreased the 
predicted velocity at the front rake a maximum of 10%. The maximum de- 
crease occurred near the augmentor wall, where the cell flow rate had the 
greatest influence. At the aft rake the 10% decrease in specified cell 
flow rate caused a maximum decrease in velocity of 7%, again near the 
augmentor wall. This decrease in specified flow rate caused the static 
pressure drop to decrease approximately 20%. These results are shown for 
the 90% power case in Figs. 28 and 29. 

The theoretical velocity profiles at the forward rake appear to be in 
good agreement with the experimental results. However, the experimental 
results showed less variations near the centerline than predicted. The 
predictions underestimated the centerline velocities at idle and 90% 
power and overestimated it at 100% power. 

The experimentally obtained aft rake velocities at idle power had 
considerable circumferential variations. The theoretical predictions, 
which are based on axisymmetric flow, appeared to be a reasonable esti- 
mation of the average behavior. In both the 90% and maximum power 
cases the computer model predicted a good velocity profile although the 
velocities were considerably greater (up to 63%) than the experimental 
values. This occurred to a greater extent at maximum power, indicating 
a possible decrease in predictive accuracy as mass flow rate was in- 
creased. Inaccurate specification of engine flow rates can significantly 
alter the cell flow rate input to the computer model. This in turn can 
Significantly affect the predicted velocity profiles as discussed above. 


Inaccurate stagnation profiles can also contribute to errors in model 
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input data as discussed above. However, these possible inaccuracies 
cannot account for the entire difference in results obtained at full 
power. 

The theoretical pressure profiles (Figs. 10, 11, and 12) in all 
three cases indicate a large pressure drop at the augmentor tube inlet 
plane. The pressure decrease was greater at the augmentor wall than 
along the engine centerline . This theoretical drop in pressure was 
greater than the experimental pressure drop, the difference increasing 
with an increase in mass flow rate. The pressure rise from the minimum 
was in better agreement with the data but also was less accurate at 
higher power settings. The predicted behavior being much greater than 
experimental results is similar to the results of Walters and Netzer 
[Ref 6 | using the t-w computer code and to the results of Stevenson and 
Netzer lRef. 1 uSing the primitive variable code. In these earlier 
Studies which compared predictions to experimental data in a one eighth 
scale test cell, good agreement in pressure profiles was obtained for 
zero engine-augmentor spacing. The discrepancy increased with increasing 
engine-augmentor spacing. One possible cause of this disagreement is the 
K-e turbulence model. It is being used in this application for the shear 
layer mixing of a practically unconfined jet surrounded by a very low 
velocity flow, both of which must then enter a confined diameter. 

Sensitivity of the computer model predictions to the magnitude of the 
effective viscosity was examined in a cursory manner by both increasing 
and decreasing the calculated values by a factor of 1.5. The computer 
model predictions were found to be fairly insensitive to these changes in 


calculated effective viscosity. For the case of a 50% increase in viscosity 
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Figs. 30 and 31 show less than 5% change in velocity at the front rake, 
15% at the aft rake, and no greater than a 6% decrease in static pressure 
drop along the augmentor wall. A 33% decrease in viscosity changed velo- 
city at the front rake less than 5%, 10% at the aft rake and decreased 
the static pressure drop along the augmentor wall no more than 9% (Figs. 
32 and 33). It appears that merely changing the magnitude of the effec- 
tive viscosity will not account for the differences in the predicted and 
measured pressure profiles. Further investigations of the two-parameter 
effective viscosity model are required. 

The sensitivity of the static pressure profiles to possible inaccura- 
cies in model inputs can account for much of the discrepancies at idle 
power but not for 90% and full power. 

Thus, it appears from both the velocity and pressure profile data 
that the predictive accuracy of the model decreases with increasing 
engine exit Mach number. 

Examination of the data from the velocimeters positioned above the 
TF-34 showed a test cell air velocity profile which increased from 
ceiling to engine. However, it has been shown by Speakman, Hayes, and 
Netzer [Ref. 5] that this condition has an insignificant effect on the 
predicted velocity and pressure profiles in the augmentor. 

The computer model appears to have the same capabilities and weak- 
nesses in predicting the flow in full scale test cells as it did for 
the subscale test cell. The subscale test cell, by design, provided 
nearly axisymmetric flow in the augmentor at all power settings and 
therefore somewhat better agreement was obtained for predicted velocity 


profiles far down the augmentor tube. Engine flow rates were also more 
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accurately known in the subscale tests. The most severe weakness of the 
model appears to be in the inability to correctly predict accurate sta- 
tic pressure profiles in regions where an expanding high velocity jet is 
captured by a fixed diameter tube. Velocity profiles and mixing rates 


appear to be reasonably predicted. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The two-dimensional, primitive-variable computer model proved to be 
quite good in predicting velocity profiles and mixing rates in the aug- 
mentor tube of a full scale jet engine test cell. The major weakness in 
this area was that the magnitude of the velocities far downstream were 
somewhat higher than actual with the discrepancy increasing with in- 
creasing engine exhaust Mach number (i.e. power setting). 

The greatest drawback to the computer model was its inability to 
accurately predict static pressure profiles in the augmentor entrance 
region where an expanding high velocity jet is captured by a fixed dia- 
meter tube. This inaccuracy also was found to increase with engine 
power level. 

If time and cost pose less restraints to future investigators it is 
recommended that the data gathering apparatus be expanded to include 
pitot-static and temperature measuring equipment at all locations where 
velocity data are desired. This would preclude the necessity for making 
the assumptions required in this investigation for stagnation tempera- 


ture profiles. 


4) 








¥3INSWONV 
90! 





suoLzeunbijuoyg 3Se, pe-di Lt “bt4 


ion seer: 
yO0O14 


Van. 
ail 


v9 G9 


ONITISS W139 


-LA INI 


42 





THERMOCOUPLE 
PROBES 


6 32.8125 


Fiad. 2 Pitot Rake Front View 
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Fig. 3 Pitot Rake Mounted in Augmentor Tube 
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Fig. 4 Pitot Rake Contruction 
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Fig. 7 Velocimeter Calibration Curves 
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Fig. 9 Static Port Labels 
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